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We carried out a feasibility study on the measurement of the branching ratio of H — > cc at a 
future e + e~ linear collider. We used topological vertex reconstruction algorithm for accumulating 
secondary vertex information and neural network for optimizing c quark selection. With an assump- 
tion of a Higgs mass of 120 GeV/c 2 we estimated statistical uncertainty of Br(H — > cc) to be 20.1% 
or 25.7% , depending on the number of vertex detector layers, at the center-of-mass energy of 250 
GeV and the integrated luminosity of 500 fb -1 . 

PACS numbers: 11.15.Ex, 11.30.Qc, 14.65. Dw, 14.80.Bn 
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I. INTRODUCTION 

In the theory of elementary particles, Higgs boson is 
introduced to elucidate the origin of particle masses. The 
Standard Model (SM) proposed the single Higgs doublet 
that gives rise to a scalar particle, while extended super- 
symmetric models hypothesized that two Higgs doublets 
give separate vacuum expectation values to the up-type 
and down-type quarks. Among the extended models be- 
yond SM, the Minimal Supersymmetric Standard Model 
(MSSM) sets the upper bound of the lightest neutral 
Higgs mass around 130 GeV/c 2 ^. Since the mass of 
SM Higgs boson is expected to be between 114.4 GeV/c 2 
and 211 GeV/c 2 by LEP experiment Q, the Higgs boson 
should be identified if it is found. Therefore a precise 
measurement of Higgs couplings to fermions at a future 
e + e _ linear collider is indispensable to unveil physics of 
the Higgs sector. In particular, a measurement of the 
branching ratio of H — ► cc would provide a unique oppor- 
tunity to study Higgs to up-type quark coupling. How- 
ever, the H — > bb is the dominant process for Higgs mass 
below 140 GeV/c 2 . 

A high efficient charm quark identification is necessary 
for the precise branching ratio study and it is possible 
only with a high performance vertex detector. Develop- 
ment of this vertex detector has been a major issue for 
the linear collider To this end, understanding detec- 
tor performance is extremely important to physics. 

In this paper, we introduce our method of measuring 
the branching ratio of H — > cc against SM e + e _ processes 
and other decay branches from Higgs. In addition, we 
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would like to see the significance of this precision study 
with four layers of Charge Coupled Device (CCD) and 
an additional inner layer of vertex detector. 



II. HIGGS SIMULATION 

A mass of SM Higgs boson is assumed to be 120 
GeV/c 2 and the center-of-mass energy of 250 GeV is se- 
lected by the highest s-channel production cross-section 
at this energy, 226 fb. It is concerned with an integrated 
luminosity of 500 fb -1 for the first few years running of 
the future e + e~ linear collider experiment. 

Events were generated using PYTHIA 5.7 [H and only 
the Higgsstrahlung process (e + e~ — > Z* — > Z°H) was 
taken into account. The branching ratios of the SM 
Higgs boson was estimated by the HDecay program 
Depending on the decay modes of Z°, the events are cat- 
egorized into 4-jet mode (Z° — ► qq and H — > qq), 2-jet 
mode [Z° — » vv and H — > qq), and charged lepton pair 
mode (Z° — > £ + £~ and H — > qq). Since the branching ra- 
tio of Higgs to cc is estimated to be very small (~3%), we 
concentrated on the 2-jet and 4-jet modes only. As back- 
ground processes from e + e~ collision, W + W~, Z°Z°, 
and qq events were considered with corresponding cross- 
sections of 15460 fb, 1250 fb, and 47300 fb, respectively 

The detector simulation was performed using a fast 
parameterized simulator [|| , which is implemented in the 
Joint Linear Collider (currently called as Global Linear 
Collider) detector Q ■ In this simulator five parameters of 
the helical track and their error matrices including non- 
diagonal elements were generated; thus the quality of the 
vertices is similar to that of a full simulation. The JLC 
vertex detector is equipped with four layers of CCD at 
the radius from 2.4 cm to 6 cm of which intrinsic spa- 
tial resolutions are 4 /jm in v<p and z directions. The 
solenoidal magnet field of the detector is 3 tesla. With 
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number of CCD layers 


Beam pipe 


1st 


2nd 


3rd 


4th 


5th 


4 


2.0 


2.4 


3.6 


4.8 


6.0 




5 


1.0 


1.2 


2.4 


3.6 


4.8 


6.0 



TABLE I: Two vertex detector parameter set. the number of 
CCD layers and their positions(in cm) from the center of the 
beam pipe, the radius of beam pipe are listed for each set. 



a vertex detector constraint, the impact parameter res- 
olution in xy plane (er^) is y/ (25/psin 2 / 3 (9) 2 + 4 2 fim, 
and the momentum resolution for charged track (Ap t /pt) 
is ^/(l x 10~ 4 p t ) 2 + 10~ 3 . This study was done with 
four/five CCD layers of vertex detector parameter, repec- 
tively, to see the influence of the vertex detector options 
on the physics result. The vertex detector parameter sets 
are shown in the Table |U 

In the event reconstruction, jets were reconstructed us- 
ing JADE clustering algorithm 0] and the vertices were 
reconstructed using the topological vertex reconstruction 
algorithm (ZVTOP program) [§]. Based on kinematic 
variables and reconstructed vertex information, an ar- 
tificial neural network (NN) [l(| is formed to identify 
H — > cc events more efficiently. 

We considered the non-cc Higgs decays (H — > 
bb, gg, WW*; Higgs background) and other processes 
from e + e~ collision (e + e~ — > Z°Z° , W + W~, qq; non- 
Higgs background) as backgrounds of H — > cc measure- 
ment. 



III. SELECTION OF H -> cc 

A. 2-jet mode 

In 2-jet mode all clustered particles are forced to make 
2 jets by adjusting the maximum value of ycut [l3l |. 
Ymax. The following selections are applied to reduce 
non-Higgs background in the sample: (1) visible energy 
between 110 GeV and 143 GeV, (2)missing transverse 
momentum between 25 GeV/c and 70 GeV/c, (3)Higgs 
mass between 105 GeV/c 2 and 125 GeV/c 2 , (4)recoiled 
Z° mass between 82 GeV/c 2 and 120 GeV/c 2 , (5)thrust 
between 0.75 and 0.99, (6)Ymax between 0.72 and 0.84, 
and (7)mass of each jet between 2 GeV/c 2 and 40 
GeV/c 2 . Here the Higgs mass is an invariant mass of 
all the observed particles and the recoiled Z° mass is 
calculated disregarding the initial state radiation. We 
also require a successful secondary vertex reconstruction 
and at least one Pt corrected invariant mass (MSPTM, 

\/M 2 vtx + Pt 2 + |Pt|) of the secondary vertex to be be- 
tween 0.1 GeV/c 2 and 7.0 GeV/c 2 , where P t is the total 
transverse momentum of the secondary tracks with re- 
spect to the flight direction of the vertex. The separation 
between signal and non-Higgs background is clearly seen 
in Fig. [TJ while Higgs background is not distinguishable 
in this aspect. 
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FIG. 1: Visible energy and Ymax distributions in 2-jet mode 
are shown. Arrows are place on the cut value. Solid: Z°H, 
dotted: W + W~, dash-dotted: Z°Z°, dashed: qq. 



After the sample selection, NN is trained to max- 
imize the selection quality. Two sets of training are 
used; one against non-Higgs backgrounds (Background 
NN-training) and the other against Higgs backgrounds ( 
Higgs NN-training). For each NN-training H — » WW* 
and e + e~ — > qq processes are not used due to the 
featureless patterns and small statistics, respectively. 

The normalized input patterns for NN-trainings were 
given from the ZVTOP program; the number of vertices, 
MSPTM, transverse momentum of the secondary ver- 
tex, decay length, invariant mass of the secondary vertex, 
number of tracks in the secondary vertex, and momen- 
tum of the secondary vertex divided by total momentum 
(corrected secondary momentum) . Distinctive input pat- 
terns for Higgs NN-training and Background NN-training 
are shown in Fig. [2] and Fig. [3J respectively. In the 
Background NN-training Ymax and jet decay angle were 
included in the input patterns especially for this 2-jet 
mode. 




FIG. 2: MSPTM and the corrected secondary momentum 
distributions are used as input patterns for H — > cc distinc- 
tion from Higgs backgrounds in 2-jet mode. Solid: H — > cc, 
dotted: bb, dashed: gg. 



Using both results of Higgs NN-training and Back- 
ground NN-training, the signal region is evaluated in two 
dimensional plane for best significance. The signal selec- 
tion area is shown in the top right corner of the Fig. [5] 
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Ymax | Decay angle 




FIG. 3: Ymax and decay angle distributions show good sepa- 
ration between cc and non-Higgs backgrounds in 2-jet mode. 
They are used as input patterns for Background NN-training. 
Solid:Z°H, dotted: dash-dotted:Z°Z°. 



| Thrust angle | | Recon. H mass. 
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FIG. 5: Thrust angle and Higgs mass distributions in 4-jet 
mode are shown. Arrows are placed on the cut value. Solid: 
Z°H, dotted: W+W~, dash-dotted: Z°Z°, dashed: qq. 
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(a) Background NN-result 
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(c) Background NN-result 



FIG. 4: 2D NN-result for H -> cc selection in 2-jet mode; 
(a)Higgs backgrounds, (b)H — > cc signal event, and (c)non- 
Higgs backgrounds. The signal region is seen at the top right. 



B. 4-jet mode 

Here we used forced 4-jet clustering for Higgs jet pair 
and Z° jet pair. The jets are identified by minimum % 2 , 
which is defined as the squared sum of the differences 
between reconstructed invariant mass of jet pair and the 
expected mass divided by each mass resolution: 
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The selection cuts were made: (1) visible energy 
greater than 210 GeV, (2) cos 9 of thrust axis(thrust an- 
de) between —0.85 and 0.85, (3) Higgs mass between 




Secondary momentum 




FIG. 6: MSPTM and the corrected secondary momentum dis- 
tributions in 4-jet mode. They are used for H — > cc distinc- 
tion from non-cc Higgs decays in Higgs NN-training. Solid: 
H — > cc, dotted: bb, dashed: gg. 



103 GeV/c 2 and 130 GeV/c 2 , (4) Z° mass between 78 
GeV/c 2 and 102 GeV/c 2 , (5) mass recoiled to the Higgs 
jet pair greater than 80 GeV/c 2 , (6) number of particles 
in each jet of Higgs pair greater than 5, and (7) Ymax 
value greater than 0.01. 

Only 2 jets which are assigned to Higgs are used for the 
vertex tagging. We also required at least one successful 
secondary vertex reconstruction with MSPTM between 
0.1 GeV/c 2 and 7.0 GeV/c 2 out of 2 jets from Higgs. The 
distributions of Higgs mass and thrust angle are shown in 
the Fig. [5j With the same procedure as the 2-jet study, 
we trained the NN against non-Higgs background and 
Higgs background separately using normalized vertex in- 
formation from the ZVTOP program. The MSPTM and 
the corrected secondary momentum distributions of the 
cc and the non-cc Higgs decays are seen in the Fig. [5] 
We did not use H — » WW* events in Higgs NN-training 
due to their featureless patterns. In the Background NN- 
training, thrust angle is added in the input pattern. As 
seen in the Fig. [7J the thrust angle shows a clear separa- 
tion between the signal and non-Higgs backgrounds. We 
use both results of Higgs NN-training and Background 
NN-training same as 2-jet mode. The 2D NN-results are 
shown in the Fig. [5] with signal box at the top right. 
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Number of CCD layers 


4 CCD layers 


5 CCD layers 


Decay mode \ Events 


2- jet 


4- jet 


2-jet 


4- jet 


H -^cc 


122.7(17.6) 


306.3(12.6) 


112.7(16.2) 


316.8(13.0) 


H -^bb 


641.8(4.2) 


2686.7(5.0) 


231.2(1.5) 


1807.2(3.4) 


H^gg 


28.1(1.8) 


217.6(3.9) 


8.9(0.6) 


143.4(2.6) 


H -> WW* 


23.6(0.8) 


226.7(2.1) 


10.1(0.3) 


178.0(1.7) 


e+e" -> W + W~ 


640(0.008) 


9790(0.130) 


330(0.004) 


8530(0.110) 


e+e~ -> Z°Z° 


100(0.016) 


1710(0.274) 


30(0.004) 


1305(0.209) 


e + e~ — > gg 


20(<0.001) 


1730(0.007) 


5(<0.001) 


1345(0.006) 


S/N 


0.0843 


0.0187 


0.1832 


0.0238 


S/VS 1 + TV 


3.09 


2.37 


4.17 


2.71 


Statistical uncertainty 


25.7 % 




20.1 % 





TABLE II: The number of c quark tagged events, efficiency (percentage in parenthesis), and significance (last three rows) 
of H — > cc measurements for two different vertex detector parameters. The statistical uncertainty (last low) is calculated 
combining significance of 2-jet and 4-jet. We used 100 fb _1 in this analysis and scaled up to 500 fb" 1 in this table. 
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FIG. 7: MSPTM and thrust angle distributions are used 
as input patterns for H — > cc distinction against non-Higgs 
backgrounds. Solid:Z°.ff, dotted: W+W", dash-dotted :Z° Z° , 
dashed :qq. 

IV. DISCUSSION 

The results of cc event selection are summarized in Ta- 
ble lU As seen in the table, roughly 120 and 310 H — > cc 
events are selected in 2-jet and 4-jet modes, respectively, 
with a reasonable significance despite its small branching 
ratio of H — > cc and large backgrounds from non-Higgs 
processes and non-cc Higgs decays. 

In the H — > cc measurement in 2-jet mode, major 
backgrounds are those from H — > bb where b is misiden- 
tified as c and W + W~ events which are reconstructed 
as 2-jet events due to imperfect detector acceptance. 
e + e~ — > evW is not considered here due to relatively 
small cross-section than other background processes. 

In the case of 4-jet mode, W + W~ and Z° Z° events are 
increased after event selection because two c jets can be 
produced in the final state of these processes. In addition 
to the c-jet contamination from backgrounds, there is 
additional ambiguity in selecting two Higgs jets out of 
four jets. Combining these effects, the significance in 
the 4-jet mode is less than 3er and worse than the 2- 
jet mode. The statistical uncertainty combining 2-jet 
and 4-jet analyses is estimated as 25.7% in the case of 
four CCD layers of the vertex detector. Similar analysis 
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(a) Background NN-result (b) Background NN-result 
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FIG. 8: 2D NN-results for H — » cc selection in the 4-jet mode; 
(a) Higgs backgrounds, (b)H — > cc signal event, and (c) non- 
Higgs backgrounds. The signal region is seen at the top right. 



with different detector configuration was report [12| for 
assumption of the same Higgs mass and luminosity at 
t/s=350 GeV and 500 GeV that the relative braching 
ratio errors be 19% and 39%, respectively. 

When an additional CCD layer of the vertex detector 
is included near the interaction point, the impact pa- 
rameter resolution for low momentum track is improved 
since the lever arm for track extrapolation to the interac- 
tion point is reduced. Thus we can reconstruct secondary 
vertices much closer to the interaction point. The bet- 
ter separation of b- and c- jets reduces backgrounds in 
the event selection, especially in 2-jet mode. Improve- 
ments in 4-jet selection is not decisive, suffering ambigu- 
ities in Higgs jet selections. Combining 2-jet and 4-jet, 
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relatively 20% improvement in background reduction is 
achieved compared to the study with four CCD layers of 
the vertex detector. 

V. SUMMARY 

In this study we focused on the H — ► cc measurements 
in 2-jet and 4-jet modes in the case of 120 GeV/c 2 Higgs 
mass at the center-of-mass energy of 250 GeV in the fu- 
ture e + e~ linear collider. In the study, the topological 
vertex finding algorithm was used for tagging the c-jet 
and the neural network was used to optimize the H — ► cc 
selection. 

With consideration of 500 fb _1 data, we obtained the 
statistical uncertainty of 25.7% for the measurement of 
H — ► cc with four CCD layers vertex detector. The sta- 
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